Ultrashort self-induced transparency plasmon solitons 
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We study the interband self-induced transmission of surface plasmon polaritons in a gold film 
surrounded by an external Kerr medium. We model the optical propagation by using a version of 
the generalized nonlinear Schrodinger equation for the field envelope coupled to Bloch equations for 
valence electrons of gold, predicting self-induced transparency of ultrashort plasmon solitons with a 
pulse duration below 10 fs. We demonstrate that the Kerr nonlinearity from the surrounding dielec- 
tric can be used to compensate for the group velocity dispersion, and that the impact of dephasing 
and decay processes can be effectively reduced by the self-induced transmission mechanism. 

PACS numbers: 42.65.Tg, 42.79.Gn, 73.20.Mf, 78.66.Bz 



Introduction - Surface Plasmon Polaritons (SPPs) 
- exponentially localized electromagnetic waves guided 
along metallic surfaces - are particularly important in 
a broad range of applications including medicine [l|, 
biosensing [2[ , imaging [3] and nanophotonics Q . A large 
variety of nonlinear mechanisms occurs in metals, e.g. 
second and third harmonic generation nonlocal pon- 
deromotive forces leading to symmetry breaking |6], and 
redshift due to interband thermo- modulation [7]. The 
subwavelength confinement of SPPs is fundamentally im- 
portant for designing nanoscaled optical circuits and for 
the enhancement of nonlinearity, which can in turn be ex- 
ploited for active all-optical control jsj, plasmon-soliton 
formation [9] and nanofocusing [To| . Inevitably, these ap- 
plications are either limited or even prevented by large 
intrinsic ohmic losses of metals. A possible strategy to 
retain the subwavelength localization of SPPs and over- 
come the loss barrier of metals at the same time is to 
embed gaining media in plasmonic devices The in- 
trinsic nonlinearity of two-level gaining media can then 
be used to obtain transverse localization of dissipative 
plasmon solitons [12]. Alternatively, losses can be lim- 
ited by using long range surface plasmon polaritons (LR- 
SPPs) in metallic films 0, 

where localization is smaller 
and absorption is reduced accordingly. In common plas- 
monic setups operating at optical frequencies, it is likely 
to use silver as metallic component in order to minimize 
losses. Indeed, while the properties of silver can be prac- 
tically understood from the free-electron plasma theory, 
the optical properties of gold and copper are greatly influ- 
enced by interband transitions, which boost absorption 
and damp optical propagation. 

In this Letter, we show the theoretical possibility to 
bypass interband absorption of LRSPPs in a thin film of 
gold surrounded by an external Kerr medium using self- 
induced transmitted 7r-pulses. Indeed, while long pulses 
(r > 100 fs) experience interband thermo-modulation 
[7[ , ultrashort pulses with time duration comparable with 
or smaller than the electron-electron collision time (r < 
T2 ~ 10 fs) do not undergo thermalization. Thus, ultra- 



short optical pulses can induce transient inversion of pop- 
ulation and experience the effect of interband polariza- 
tion. If one neglects non-conservative dephasing and re- 
combination processes, perfect self-induced transparency 
(SIT) can be achieved with 7r-pulses [14j. Moreover, for 
some specific parameters, SIT pulses can coexist with 
temporal plasmon solitons, where the group velocity dis- 
persion (GVD) is compensated for by the Kerr nonlinear- 
ity, analogously to what happens in erbium-doped fibers 
[15j. As a consequence of non-conservative processes, 
ideal transparency can not be achieved but enhanced self- 
induced transmission occurs, as we show here for the first 
time. A sketch of the setup considered in our calculations 
is depicted in Fig. [TJ 

Background - The temporal dynamics of electrons 
in the presence of an external electromagnetic field is 
determined by the Schrodinger equation with unper- 
turbed Hamiltonian T-Lq and dipole interaction Hamilto- 
nian Hi(t). Since the electron momentum is conserved 
in the dipole interaction, it is possible to split the di- 
electric polarization of gold in intraband (I) and inter- 
band (II) contributions. For intraband transitions, due 
to the high occupation of states in the conduction band, 
the periodic potential of the lattice is screened by the 




FIG. 1: Sketch of the representative structure analyzed in 
this work. Gold film of width a = 20 nm surrounded by an 
external Kerr medium. The structure is assumed infinitely 
extended in the z,y directions. 
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FIG. 2: (a) Real and (b) reciprocal lattice structures of gold. 
(c,d) Band structure around the (c) L and (d) X points of 
the reciprocal lattice. The upper blue and lower red curves 
are the conduction and valence bands, respectively. The full 
black lines represent the Fermi level E — Ef = 0. Note that 
the right and left sides of the x-axis correspond to the lattice 
vector components k±,k// rescaled to kr = 13.34 nm -1 . 



quasi-free conduction electrons. The unperturbed Hamil- 
tonian can then be approximated by the kinetic part 
Ho ~ p 2 /(2m e ) and thus the dielectric polarization is 
Pj = eoX/E, where eo is the dielectric permittivity of 
vacuum, \i = -Up/[u(w + ir )L = 13 - 8 x 1q3 THz is 
the plasma frequency of gold, T = 100 THz is the damp- 
ing rate due to electron-electron collisions, E = Eoe _lu;£ 
is the electric field and uj is its angular frequency. 

Interband transitions and Block equations - The real 
and reciprocal lattice structures of gold are face cen- 
tered cubic (f.c.c.) and body centered cubic (b.c.c), re- 
spectively, as shown in Figs. [2ja,b). Experimental re- 
sults clearly indicate that interband absorption is reso- 
nant around the X, L points of the first Brillouin zone 
[16| [see Fig. Efb)]. Around these band-edge points the 
Fermi surface is cylindrically simmetric, the lattice vec- 
tor k can be split into in-plane (k±) and out-of-plane 
(k//) components and the valence E v (k) and conduction 
£? c (k) bands can be approximated by elliptic and hyper- 
bolic paraboloids [7]. In the following calculations, we 
approximate the conduction band to be completely filled 
up to the Fermi level E < Ep and completely empty 
above it E > Ep. The band structure around the L,X 
points of the reciprocal space is displayed in Figs. [2fc,d). 
The upper blue and lower red curves are the conduction 
and valence bands, while the full black lines indicate the 
Fermi level E = Ep = 0. Yellow circles represent occu- 
pied states, while full black (dashed cyan) arrows indicate 
direct interband transitions allowed (forbidden) by the 
Pauli exclusion principle. The modeling of electron inter- 



band dynamics is developed as a first order perturbative 
solution of the Schrodinger equation in the Heisenberg 
picture, calculating the time-dependent inversion of pop- 
ulation and the complex polarization from Bloch valence 
|i7,k) and conduction |c, k) states. In our analysis, we 
neglect the k-dependence of the dipole matrix element, 
which is approximated by d CiV (k) = (c, k|er|u,k) ~ d CiV . 
Thus, the interaction Hamiltonian is explicitly given by 

W/,kW = -S(t){dc,t;|c,k)(v,k|+^hk)(c,k|}, (1) 

where we have used d* v = d V}C . From the equation of mo- 
tion for the density matrix p\z(t) = — (i / h)[H\c(i) , Pk(i)], 
one gets the Bloch equations for the population difference 
= p c ,c(k, t) — p ViV (k,t) and for the complex dipole 
P\a — d v ^ c p c ^ v (\<.jt): 

w k = -^r(l + w k ) + ^E{t)lm{P k }, (2) 

1 1 fl 

Pk = -^Pk-i^Pk- ^E(t)\d v , c \ 2 w k , (3) 

where ^ k = [E c (k)-E v (k)]/h and 7\ = 500 fs, T 2 = 10 fs 
are the decay and dephasing times of gold. The complex 
interband polarization P/j is parallel to the electric field 
E and its amplitude is Puit) = (27r 3 ) -1 J q Pkd 3 k, where 
the k-integral is taken in the volume ft of the reciprocal 
space where the interband transitions are not forbidden 
by the Pauli exclusion principle. 

Ultrashort SPP pulses in a gold film - We model 
optical propagation of SPPs starting from the time- 
dependent Maxwell equations for the electric E and mag- 
netic H fields. Since the dielectric susceptibility of gold 
can be separated in intraband and interband contribu- 
tions, it is possible to construct the linear modes with 
the intraband polarization P/ and study their nonlinear 
evolution due to the interband counterpart P/j. The 
space-dependent linear susceptibility profile of a thin 
gold film of width a surrounded by an external dielectric 
medium is €l(x,oj) = €d(oj)0(\x\— a/2)+ei(uj)9(a/2 — \x\), 
where Q(x) is the Heaviside step function, €i(uj) = 1 + 
Xi(w) is the intraband dielectric constant of gold and 
€d(uo) is the external dielectric constant. The space- 
dependent nonlinear polarization is Pnl = ^n0{a/2 — 
\x\) + (eoX3/2) [|E| 2 E + E 2 E*/2] 6{\x\ - a/2), where X s 
is the Kerr susceptibility of the external medium (xs ~ 
2.25 x 10 -22 m 2 /V 2 for silica glass) and P/j is the com- 
plex interband polarization of gold (implicitly nonlinear) 
calculated from the Bloch equations. In order to derive 
an effective propagation equation for the field envelope, 
we proceed in a perturbative manner. First, we calculate 
the linear modes of the plasmonic structure neglecting 
the nonresonant losses (e") and the nonlinear polariza- 
tion (Pnl)- Then, we perturbatively derive a propaga- 
tion equation for the linear field envelope accounting for 
nonresonant losses, Kerr nonlinear ity and interband po- 
larization. As the structure is assumed to be infinitely 
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FIG. 3: (a) Longitudinal (e z ) and transverse (ie x ) SPP mode 
profiles, (b) SIT plasmon soliton peak power density l^ol 2 as 
a function of its time duration r. In both figures the thickness 
of the gold film is a = 20 nm, the external medium is silica 
glass and the optical wavelength is A = 623 nm. 
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extended along the ^-direction, the modal electric E and 
magnetic H fields are E = y/2iioc/ S z i/je(x)e il3oZ ~ iu!ot and 
H = y/2/bi cS z )iljh(x)e i/3oZ - iuJot . In the Ansatz above, 
the fields are normalized to S z = Re[e x h*] • zdx, 
/io,c are the vacuum magnetic permeability and speed 
of light, cjo,A) are the modal angular frequency and 
wavenumber, e = e x x + e z z and h = (c/uo)(Poe x + 
ide z /dx)y are the dimensionless electric and magnetic 
profiles. The modal normalization constants are chosen 
in such a way that |^| 2 is the power density (measured 
in W/m) carried in the z-direction. 

The electric field components e Xl e z and the propa- 
gation constant /3q can be calculated by substituting 
the Ansatz in the time-dependent Maxwell equations 
and by imposing the boundary conditions (BCs) for 
the continuity of the longitudinal electric field (e z ) and 
transverse displacement vector (e' L e x ) at the interfaces 
x = d=a/2. By doing so, one gets the dispersion re- 
lation of LRSPPs tanh(g m a/2) = —qdc'i/(qmCd)i where 
Vim = Pi ~ ^e^j/c 2 . The LRSPP longitudinal (e z ) 
and transverse (ie^) field profiles are shown in Fig. [3^a). 
In the thin film limit q m a <C 1, the transverse electric 
field profile (e x ) is approximately constant and the lon- 
gitudinal electric field profile (e z ) is approximatley lin- 
ear. While the transverse field component is never null 
inside the metal and stays roughly constant, the longi- 
tudinal field component vanishes exactly at the center 
of the gold film. For this reason, the major role in the 
interband dynamics is played by the transverse field com- 
ponent that within the gold film can be approximated by 
e x ~ — 2i/?o/((7m a )- I n t ne following step of the pertur- 
bative analysis, the mode envelope ip(z,t) is allowed to 
evolve on a slower scale compared to the fast oscillations 
e i(3oz-iu, t : 1^ ^ |^| ^ |^ o ^ Developing 

a multiscale expansion and applying the Fredholm al- 



FIG. 4: Conservative temporal dynamics of hyperbolic se- 
cant pulses in absence of intraband loss, dephasing and de- 
cay. (a,b,c) Temporal evolution of (a,b) tt and (c) 1.67T pulses 
of duration r — 5 fs in a thin film of gold of width a — 20 
nm surrounded (a) by silica glass (7 = 3.75 x 10 -7 W _1 ) and 
(b,c) by a defocusing medium that satisfies the Nakazawa con- 
dition (7 = -2.99 x 10~ 5 V^ _1 ). (d) Evolution of pulse area 
for several input hyperbolic secant pulses with different in- 
put amplitude. Note that while the area of the SIT plasmon 
soliton remains constantly 7r, for other input pulses the pulse 
area does not remain constant. 



+ 7 |V| 2 V' + itp / p k d 3 k = 0, 



ternative theorem with nonlinearly reinforced BCs [12] 
it is possible to prove that one obtains the following 
generalized nonlinear Schrodinger Equation (GNLSE) for 



(4) 



where £>k = P k e _1 ^° 2+lu;ot , v g = du/d/3 is the group ve- 
locity, /?2 = d 2 f3/duj 2 is the dispersion coefficient, a is 
the absorption coefficient due to intraband losses, 7 is 
the Kerr nonlinear coefficient and ft p is the interband 
coupling coefficient. The formulation of Eq. dU is the 
first major result of this Letter. 

Self-induced transmission of SPPs - The complex in- 
terband polarization £>k couples the GNLSE with Eqs. 
(03]), which under the rotating wave approximation 
(RWA) become 

^ k = C[# k + ^*Pk] - (1 + ufc)/Ti, (5) 
p k = iA k Pk - (C/2)\d V:C \ 2 i/;wk -pk/T 2 , (6) 

where Ak = ujo — is the k-dependent detuning and 
C = 2/3o^ s /2j~Loc/ (q 2 n ah\fSl~) is the interband coupling con- 
stant. For pulses with time duration r much smaller 
than decay (Ti) and dephasing (T2) times, the non- 
conservative terms of Bloch equations can be disregarded. 
In this limit, neglecting also the intraband loss (a), it 
is possible to derive the expression for the family of 
self-induced transparency (SIT) solitons [14j: ip(z,t) = 
^ sech[(t-z/V)/r]e iA ^, where ^0 = 2 (C|4, c |r) _1 , V 
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is the SIT soliton velocity and A/3 is its phase shift. The 
dependence of the SIT plasmon soliton power density 
( | ipo | 2 ) on its time duration (r) is depicted in Fig. [3fb). 
Assuming that the lateral confinement in the unbound 
^-direction is of the order of L y w 10 /im, one gets that 
SIT plasmon solitons can be excited with peak powers 
of the order of P = \ipo\ 2 L y ~ 10 kW. We emphasize 
that the SIT plasmon soliton exists only for a particu- 
lar coincidence of geometrical parameters such that the 
GVD is exactly compensated for by the Kerr nonlinear- 
ity faC 2 \d V:C \ 2 + 47 = 0. The interband absorption is 
resonant at A = 400 nm and has a characteristic spectral 
width of AA ~ 100 nm Q. Thus, for A > 500 nm the 
band dispersion does not play a crucial role, it is possi- 
ble to approximate Ak ~ A and the interband dynamics 
is simplified to a two-level system with effective detun- 
ing A. In this approximation, we have solved numeri- 
cally Eqs. (|4|5|6|) by using the split-step Fourier method 
with fourth-order Runge-Kutta algorithms for the spa- 
tial propagation and the Adams-Bashforth method for 
the temporal integration of Bloch equations. In the sim- 
ulations, we considered a thin film of gold with thick- 
ness a = 20 nm surrounded by silica glass at the optical 
wavelength A = 623 nm, characterized by group veloc- 
ity v g = 194 km/s and dispersion coefficient fa = 0.92 
/im -1 fs 2 . Results of numerical simulations in the conser- 
vative limit are reported in the panel of Fig. HI In Figs. 
2](a,b,c) we show the temporal dynamics of hyperbolic se- 
cant pulses of duration r = 5 fs with area (a,b) tt and (c) 
1.67T in a thin film of gold of width a = 20 nm surrounded 
(a) by silica glass (7 = 3.75 x 10" 7 W" 1 ) and (b,c) by 
a defocusing medium that satisfies the soliton condition 
existence faC 2 \d VyC \ 2 + 47 = 0, with 7 = -2.99 x 10~ 5 
W _1 . In Fig. Hfd), the evolution of the pulse area is 
depicted for several input hyperbolic secant pulses with 
different amplitude. Note that, when the soliton exis- 
tence condition is satisfied, the propagation of tt pulses 
does not affect their shape, while pulses with different 
areas are either absorbed or undergo an oscillatory dy- 
namics. The discovery of such ultrashort SIT solitons in 
gold is the second major result of this Letter. 

Results of numerical simulations in realistic non- 
conservative conditions are reported in the panel of Fig. 
[U In Figs. [H[a,b,c) we show the temporal dynamics of 
hyperbolic secant pulses of duration r = 5 fs with tt area 
in a thin film of gold surrounded of width a = 20 nm 
(a) by silica glass (7 = 3.75 x 10 _7 W _1 ) and (b,c) by 
a defocusing medium that satisfies the soliton condition 
existence faC 2 \d V)C \ 2 + 4 7 = 0: 7 = -2.99 x 10~ 5 W" 1 . 
In Fig. Sfd), the evolution of a long pulse of duration 
r = 100 fs is depicted. In this limit thermalization oc- 
curs and the optical propagation is modeled in terms of 
an effective dielectric constant. Note that, even if ideal 
transparency is not achieved, ultrashort pulses of dura- 
tion r ~ 10 fs propagate over much longer distance with 
respect to long pulses of duration r ~ 100 fs. 
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FIG. 5: Non- conservative temporal dynamics of hyperbolic 
secant pulses. Temporal evolution of tt pulses of durations 
(a,b) r = 5 fs, (c) r = 20 fs and (d) r = 100 fs in a thin 
film of gold of width a = 20 nm surrounded (a) by silica glass 
(7 = 3.75 x 10 -7 W _1 ) and (b,c,d) by a defocusing medium 
that satisfies the soliton existence condition (7 = — 2.99x 10 -5 

w- 1 ). 



Conclusions - We have investigated analytically and 
numerically the interband self-induced transmission of 
surface plasmon polaritons in a gold film, demonstrat- 
ing the possibility to bypass the interband losses of gold 
by using ultrashort 7r-pulses. This regime is realistically 
accessible to experiments with pulse peak powers of the 
order of P ~ 10 kW and pulse time duration r < 10 
fs. These theoretical findings highlight the novel possi- 
bilities offered by the use of ultrashort pulses in plas- 
monic systems and can stimulate further developments 
and strategies to reduce the impact of losses in metals, 
as well as opening up new directions for the exploitation 
of nonlinear effects in metals. 



* Electronic address: andrea.marini@mpl.mpg.de 
[1] A. M. Gobin, M. H. Lee, N. J. Halas, W. D. James, R. 

A. Drezek and J. L. West, Nano Lett. 7, 1929 (2007). 
[2] J. N. Anker, W. P. Hall, O. Lyandres, N. C. Shah, J. 

Zhao and R. P. Van Duyne, Nat. Mat. 7, 442 (2008). 
[3] S. Kawata, Y. Inouye and P. Verma, Nat. Phot. 3, 388 

(2009) . 

[4] D. K. Gramotnev and S. I. Bozhevolnyi, Nat. Phot. 4, 83 

(2010) . 

[5] P. Ginzburg, A. V. Krasavin and A. V. Zayats, New J. 
Phys. 15, 013031 (2013). 

[6] A. R. Davoyan, I. V. Shadrivov and Y. S. Kivshar, Opt. 
Lett. 36, 930 (2011). 

[7] A. Marini, M. Conforti, G. Delia Valle, H. W. Lee, Tr. X. 
Tran, W. Chang, M. A. Schmidt, S. Longhi, P. St. J. Rus- 
sell and F. Biancalana, New J. Phys. 15, 013033 (2013). 

[8] M. Kauranen and A. V. Zayats, Nat. Phot. 6, 737 (2012). 



5 



[9] E. Feigenbaum and M. Orenstein, Opt. Lett. 32, 674 
(2007). 

[10] A. R. Davoyan, I. V. Shadrivov, A. A. Zharov, D. K. 
Gramotnev and Y. S. Kivshar, Phys. Rev. Lett. 105, 
116804 (2010). 

[11] M. C. Gather, K. Meerholz, N. Danz and K. Leosson, 

Nat. Phot. 4, 457 (2010). 
[12] A. Marini and D. V. Skryabin, Phys. Rev. A 81, 033850 



(2010). 

[13] P. Berini, Adv. in Opt. and Phot. 1, 484 (2009). 

[14] S. L. McCall and E. L. Hahn, Phys. Rev. 183, 457 (1969). 

[15] N. Nakazawa, E. Yamada and H. Kubota, Phys. Rev. 

Lett. 66, 2625 (1991). 
[16] N. E. Christensen and B. O. Seraphin, Phys. Rev. B 4, 

3321 (1971). 



